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SUMMARY 


r 

A  computational  algorithm  utilizing  the  free  vibration  modes  of  a 
structure  is  presented  for  calculating  receptances.  The  usual  eigen- 
system  computed  for  large  structural  models  is  incomplete;  hence  the 
receptances  are  approximate.  The  formulas  developed  here  increased 
receptance  accuracy  compared  to  classical  spectral  representations. 
Receptances  are  used  extensively  in  eigensolution  reanalysis,  design 
and  synthesis  and  also  for  forced  harmonic  response  studies.  In  these 
areas  receptance  approaches  offer  a  popular  alternative  to  Rayleigh 
Ritz  subspace  type  methods.  Structural  models  represented  by  non- 
syrametric  mass,  damping  and  stiffness  matrices,  which  occur  frequently 
in  rotating  structures,  may  be  treated  using  the  receptance  formulas 
presented.  The  receptance  matrix  derived  is  applicable  to  general, 

2nd  degree,  square  lambda  matrices.  This  generalizes  the  receptance 
matrix  commonly  associated  with  matrix  pencils. 
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jounce, 
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INTRODUCTION 


In  the  terminology  of  vibration  theory  the  receptance,  or 
dynamic  flexibility  influence  coefficient,  between  degrees  of  freedom 


i  and  j  is  defined  as 

G_(u>)  =  complex  amplitude  of  the  steady  state  displacement 

at  dof  i  due  to  a  unit  sinusoidal  force  of  frequency 
u>  applied  at  dof  j  . 

=  {(-U)2  M  +  iu>  C  +  (1) 

Receptances  have  been  extensively  used  in  dynamic  reanalysis,  that  is 
in  the  efficient  solution  of  the  matrix  eigenvalue  problem  for  a  struc¬ 
ture  undergoing  physical  modifications.  Receptances  are  also  used  in 
optimized  design  and  structural  damage  studies  as  well  as  in  the  syn¬ 
thesis  of  structures  possesing  prescribed  eigenvalues.  In  the  latter 
case  this  procedure  is  referred  to  as  eigensolution  design  or  "pole 
placement."  Modal  synthesis,  i.e.,  the  computation  of  eigensolutions 
of  a  segmented  system  from  the  eigensolutions  of  its  component  substruc¬ 
tures,  has  also  been  presented  from  a  receptance  standpoint.  The  recep¬ 
tance  methods  can  in  some  cases  provide  much  better  computational 
efficiency  and  accuracy  than  a  standard  Rayleigh  Ritz  subspace  approach. 
In  reanalysis  this  is  particularly  evident  if  the  subspace  dimension 
of  the  Rayleigh  Ritz  method  is  much  larger  than  the  number  of  degrees  of 
freedom  directly  affected  by  the  physical  modification. 

Bishop  et.  al.  systematically  utilized  receptances  in  vibration 
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theory  and  experimentation.  ’  ’  To  calculate  individual  elements  of  a 
receptance  matrix  for  an  undamped  system,  these  authors  utilized  its 


spectral  representation 


This  relationship  permits  any  element  of  the  receptance  matrix  to 
be  computed  at  various  forcing  frequencies  (to) ,  without  repeatedly  in¬ 
verting  large  matrices.  Since  the  modal  characteristics  calculated  for 


complex  structural  models  are  usually  incomplete,  the  summation  in 
equation  (2)  extends  only  over  the  L  lowest  modes.  Typically, 

0.1  <  L/N  <  0.3  (3) 

Hence  the  spectral  representation  of  the  receptance  matrix  is  aporoximate. 

Receptance  investigations  have  sought  to  improve  the  representation 

of  equation  (2)  for  L/N  values  in  the  range  of  equation  (3).  Hirai  and 
4 

YosLimura  developed  an  alternative  receptance  formula  useful  in  eigen- 
solution  reanalysis  of  locally  modified,  structural  systems.  Their 
"hybrid"  representation  is 


G(u>)  =  {IfJ+U)2K  1M+uj4(K'1M)2+. 


•+(w2)R"1(K~1M)R"1}k' 


+  uj' 


2R 


N  *.  £ 

y  _ J J - 

•  i  —  y-»2R,-r,2  2\ 

1=1  m.Q.  (Q.-ur) 
J  J  J  J 


(4) 


for  R  »  1,2, .. . 

They  demonstrated  that  the  accuracy  of  G  (u>)  improves  as  R  increases, 
for  a  fixed  value  of  L.  Leung^  utilizes  receptances  in  forced  harmonic 
response  analyses  and  eigensolution  reanalysis  of  structural  systems 
subjected  to  modified  boundary  conditions.  His  improvement  to  the 
receptance-spectral  representation  of  equation  (1)  resembles  the  R  =  2 


form  of  equation  (4) . 

Irabert^  separates  equation  (2)  into  two  series 


ha  ^2  +  £ 


G .  .  (u>)  =  2  — J--  -  +  I  -  z 

1J  £=1  ma^Qa  '  uj2-)  £=L+1 


Assuming  that  u)  <<  equation  (5)  becomes 


L  *i£  *i£  -1  L 

G.  .  (u))=  2  .-■J  -  Mj!  )..  •  2  -±=— ^ 

1J  £=1  m£(n|-u)2)  1J  a- 1  m£Q2 


-1  ^  ^i£  ^i£ 

(K  1 ) .  .  =  2  L» 

1J  £=l  “£Q£ 


It  can  be  shown  that  equation  (6)  and  the  P  =  1  form  of  equation  (4)  are 
equivalent.  The  last  two  terms  in  equation  (6)  are  referred  to  as  the 
residual  flexibility,  and  take  into  account  the  static  effects  of  the 
missing  higher  order  modes.  Childs7,  in  a  study  of  transient  modal 
rotordynamic  models,  has  recently  applied  the  residual  flexibility 
approach. 

The  present  work  deals  with  the  more  general  receptance  formed  as 
the  second  degree,  lambda  matrix  inverse 

F(\)  =  (\2M+AC+K)_1  NXN  (8) 

where  M,  K,  and  C  are  square,  generally  nonsvmmetric  matrices.  The 
scalar  \  is  in  general  a  member  of  the  complex  number  field.  The  matrix 
F  will  be  referred  to  as  a  generalized  receptance  matrix.  The  generalized 
receptance  is  useful  in  vibration  theory  while  examining  rotating  struc- 

Q 

tures  ,  structures  with  rotating  elements,  and  structures  with  discrete 
dampers.  This  last  category  is  gaining  increasing  importance  due  to  a 


growing  interest  in  active  vibration  control  systems 


MATHEMATICAL  FORMULATION 


We  begin  the  study  of  equation  (8)  by  considering  the  quadratic 


eigenvalue  problem 


(a2  M  +  a.  C  +  K)  A.  =  0 

1  -  i  -  -  -i 


(NX1) 


and  the  corresponding  left  eigenvector  problem 


(tt2  MT  +  Of.  CT  ♦  KT)  6.  =0 
1  -  1-  -  “l  “ 


(NX1) 


These  equations  are  expressed  in  first  order  form  asJ 


<“k  ii  *  8iUk  =  ?  ■  <«k  si  *  $I>rk  =  ° 


M  '  0 


A  - a - ,  B  -  - + - 

0  -K  1  K  I  0 


C  \  K 


or  alternatively  as 


h  *  52Kk  =  o  .  (\  -  o 


0  J  m" 

A  = - + - 

M  J  C 


-M  1  0 

i  - 

0  i  K 


The  modal  matrices  associated  with  equations  (9)  -  (14)  are 


9  -  liji  A2i  ...  !A2N] 


!A2n]  ,  »  «  l«,  I  «2  ; 


r  -  rr  !  r  !  T  1  =  -  — 

L  ,  L2:  •  •  ■  ;L2nj  _ i. 

_  9 


r  *  h  i  U2:  iW 


_  J  V  diag  (or.) 


f 


Assume  that  each  eigenvalue  Of^  is  distinct  with  an  eigenspace  of  dimen- 
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sion  one.  The  oiorthogonality  conditions  satisfied  by  these  eigen¬ 
vectors  then  are 


(v^A,  T) . .  =6..  a..,  (v^  B.  r)..  =  -6..  a  a. 

-1  -  n  ii  li’  -1  -1  - H  lj  li  i 


(18) 


(\T  A,,  H  .  .  =  6.  .  a~ .  ,  (vT  B„  n . .  =  -6. .  a  a. 
-2  -  ij  ij  2i  ’  A  -2  -  ij  ij  2i  i 


(19) 


a..  =  6.(2a.M  +  C)A.  ,  a  =  a. a„ . 

2i  -ii-  —-ill  i  2i 


a^O 


(20) 


Define  the  matrix  pencil 


P  =  A  A,  +  B,  ,  \  f  cr 


(21) 


The  inverse  of  this  pencil  is  obtained  using  equations  ( 1 6 )  -  (20) 


Substitute  equations  (14),  (21),  (22)  into  the  identity 


P  p'1  :  I 

-  -  -2N 


(23) 


This  yields  the  conditions 


-M  Q  diag  f  °i  \  +  A  M  Q  diag  f  1  \ 

"  ""  la2i(^i  V 


V  =  0  (24) 


AM  3  diag  f  ai  +  (AC+K)Q  diag  _ 1 _ _ -  1^ 


.32iiX-aiV 


La2i(V“il 


(25) 


f 


Multiplying  equation  (24)  by  A  and  adding  the  results  to  (25)  provides 
the  spectral  representation  of  the  generalized  receptance  matrix 


2N  A. 6. 

F(A)  =  I - (26) 

j=l  a2j(A-a^) 

When  the  structure  is  constrained  and  unable  to  undergo  a  rigid  body 
displacement  in  any  vibratory  coordinate  direction,  all  2>J  eigenvalues 
(o^)  are  non-zero.  In  the  limit  as  A  approaches  zero  equation  (8)  and 


(26)  imply 


.  2N  A.  6 
K_I  =  -  X 

j=l  »2J  a 


Consider  the  following  finite  expansion 

1  1  *  A  A2  A3 

A-a .  a .  ,  "  "  “  " 

J  J  a2  a? 

3  )  1 

_  A*'1  ,  AR 

aR  c*R(,\-of . ) 
J  J  J 


R  =  1 ,  2 ,  . . . ,  R  finite 


Substitute  equation  (28)  into  (26) 


2N  A.  ST  2N  A.  6T 


2N  A.  6 

F(A)  =  -  I  J-  -  A  i  -J - A  -  A2  I 

j  =  1  2J  j  j  =  l  a2i  j  =  1  a?4  a-3 

.1  2J  1 


2N  A.  ST 


-A3  I 


j=l  32j  a] 


B  ,  2N  A.  6T  _  2N  A.  5T 

..  -  AR_1  I  +  Z  - - 

j=l  a.  .  qR  j  =  l  a0.  aR(A-cf.) 

2j  J  J  2j  j  j 


It  follows  from  equation  (27)  that  the  first  term  in  this  series  equals 


the  inverse  of  the  stiffness  matrix.  Likewise,  the  first  R  terms  of 


(29)  may  be  expressed  in  terms  of  the  M,  K  and  C  property  matrices. 
These  expressions  are  obtained  by  equating  the  spectral  representation 
for 


Bj1  (At  B^)R'2 


R>2 


(30) 


to  its  representation  as  defined  by  equations  (12).  The  spectral  repre¬ 
sentations  for  and  are  derived  using  equations  (18)  -  (20). 

The  spectral  representations  for  and  are  derived  using  equations 
(18)  -  (20). 


•V 


=  f  diag 


a .  a  . 

J  2j, 


a2  A . 

.T  1 
o  .  i  a  . 

A. 

„T  “1 
o . 

J  -J 

-J  |  J 

-J 

-J 

— 

—  -+  — 

— 

— 

2N 

a .  A . 

«.T  1 

0  .  1 

A. 

6T 

-  y 

-  J  “J 

-J  1 

—  ] 

-J  _ 

jh 

a .  a 

J  2j 

(31) 


Si 


r  diag  j  1 


a2  a  . 

"  J  -J' 


2N 

-  I 

j  =  l 


a. 

T  1 

•5  i  a.  A. 

j  -j 

-J 

J  -J 

“J 

— 

- 4 _ 

_  — 

a .  A . 

T  1 

6l  i  A. 

6T 

L  LZJ. 

“J 

-J  J 

n2 

*J  "21 


(32) 


In  deriving  equivalent  expression  in  terms  of  M,  jK,  and  C  for  the 
summations  in  equation  (29),  only  the  "2,2"  submatrix  of  equation  (30) 
need  be  examined.  The  hybrid  spectral-property  matrix  representations 
for  the  generalized  receptance  matrix  for  R  equal  to  1  through  4  are 


R  =  0 
F(A)  = 


2N 

I 

j  =  l 


A.  6T 

-J  'J _ 

a2j(A-a.) 


J 


ii 


(33) 


R  =  1,  «  .  t  0 


IMPROVING  THE  COMPUTATIONAL  ACCURACY 


The  integer  R  in  equations  (29)  and  (33)  -  (37)  can  be  referred  to 
as  a  "higher  mode  factor."  To  investigate  why  increasing  R  reduces  the 
importance  of  the  higher  mode  contributions  to  the  summations  in  equations 
(33)  -  (37)  consider  the  ratio 


T 

coefficient  of  A  6 
-n  -n 

R 

a91(A  -  ffj) 

T 

coefficient  of  5^ 

% 

a9  (A  -  a  ) 
2n  n 

Arrange  the  eigensolutions  according  to 


“l  -  a2 


<  a. 


w  < 

2N-1  - 


2N 


(39  > 


where  the  equality  holds  only  if  +  ^  -  ay  This  implies  that 


<  1  ,  3  <  n  <  2N 


(40) 


Consequently  equation  (38)  shows  that  as  R  increases,  r_^  decreases  for  n 
greater  then  three  (3)  .  Hence  the  higher  mode  contributions  to  the 
summations  in  equations  (33)  -  (37)  are  increasingly  less  important  for 
larger  values  of  R.  Since  these  modes  are  not  usually  calculated  in  the 
analysis  of  large  structure  models  the  accuracy  of  the  elements  of  F(A) 
will  improve  as  R  increases. 


Consider  the  contribution  of  the  mth  mode  to  the  rec  ptance 
matrices  in  equations  (33)  -  (37),  i.e., 

A  5T 

—to  — m  _ 

a9  (A-ot  ) 
dm  ra 


R  =  1 


A 


A 

ra 


6 


T 

m 


a  a„  (A-d  ) 
ra  2m  ra 


X 

a 

ra 


G 

— m 


R  =  2 


i  \  6 1  ■> 

A1  _  =  (1-)"  «■ 

a2  a2m(X_am)  am  ^ 
m 


Clearly  the  magnitude  of  the  contribution  from  the  mth  mode  will  de¬ 
crease  as  R  increases  if 


—  <  1 

ot  I 

m  1 


This  shows  that  if  only  q  ei^ensolut ions  are  utilized  in  equations 
(33)  -  (37)  then 


~  I  <  i 
■  \  ! 

is  a  necessary  condition  tor  the  R  =  1,2,3...  receptance  formulas  to  produce 
more  accurate  results  than  the  R  =  0  formula.  Therefor^,  the  (R  <_  1)  receptance 
formulas  for  j[(X)  should  only  be  used  if  jXj  <  |a  ]. 

APPLICATION  TO  SPECIAL  SYSTEM  TYPES 


The  generalized  receptance  formulas  given  by  equations  (29),  and 
(33)  -  (37)  can  be  simplified  for  several  common  types  of  vibrational 
systems.  Gyroscopic  systems  satisfy 


M  >  0,  K  >  0,  C  =  -C  (41) 

14,15 

These  frequently  occur  in  models  of  rotating  structures  ’  .  For  this 
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system  type 


6.  =  fty  Re  (a.)  =  0,  Refa^)  =  0 


r  "-V .  *v 


•  litnteigiiri  m 


The  eigenvectors  of  general,  symmetric  systems 


Consider  a  system  model  where 


(kj  k2  k3  k4)  =  C j t  j2  J3  J4)  =  (6  7  8  9) 
The  R  =  1  formula  from  equation  (34)  implies 


«-H  VO 

**? 

I _ 

K'1 

-67 

K_1 

^68 

-69 

466 

A67 

A68 

A69 

f(A)  = 

K'1 

-76 

K'1 

—77 

K'1 

-78 

-79 

+\ 

K"1 

-86 

-87 

K"1 

-88 

k"1 

-89 

A86 

A87 

A88 

A89 

K_1 

-96 

-97 

K'1 

-98 

K'1 

-99 

A96 

A97 

A98 

A99 

Ai  =  l  _ 

ij  £-1  a2£  a£(\  -  cr£) 

The  R  =  2  formula  from  equation  (35)  becomes 


-66 

K'1 

-67 

K_1 

-68 

f(\)  = 

-76 

K'1 

-77 

K"1 

-78 

K*1 

-79 

K-1 

-86 

K'1 

-87 

K"1 

-88 

-89 

K-1 

-96 

k"1 

-98 

K99 

A66 

4 

*?7 

4 

‘■h 

A87 

A88 

A8' 

A96 

A97 

A98 

4 

where 


1 


'V-  1-1  -  1 

and  K  and  K  are  the  nth  row  and  column  of  K  ,  respectively, 
-n  — n  ~ 

Similar  formulas  for  R  greater  than  2  are  also  easily  derived.  Note 
that  the  property  matrix  -  submatrices  in  equations  (49)  and  (50)  need 
be  evaluated  and  stored  only  once.  These  submatrices  may  then  be  re¬ 
trieved  to  evaluate  f(\)  for  any  value  of  K. 

APPLICATION  OF  THE  RECEPTANCE  FORMULA 


A  useful  application  of  receptance  formulas  is  in  dynamic  reanalysis. 
Suppose  that  modifications  to  a  structure  result  in  the  formation  of  the 
new  mass,  damping,  and  stiffness  matrices 

M  =  M  +  AM  (N  x  N)  (51) 

C  =  C  +  AC  (52) 

K  =  K+AK  (53) 

The  quadratic  matrix  eigenvalue  problem  for  the  modified  system  becomes 


( X  (M  +  AM)  +  X  (C  +  AC)  +  (K  +  AK)  )'i>  =  0 

n  —  —  n  —  --  —  —  -n  — 


(54) 


From  equation  (54)  define  the  modification  matrix 

S(X)  =  A2  AM  +  AC  +  AK  (55) 

Suppose  that  the  modification  involves  relatively  few  elements  in  the 
property  matrices,  in  other  words,  AM  etc.,  is  very  sparse.  The  sparse¬ 
ness  of  the  modification  is  quantified  by  considering  the  non-null  rows 
(or  columns)  of  S^(A).  These  row  numbers  are  entered  in  the  set 

1  =  (j i  j 2  ’ • •  Jp)T  P  x  1  (56) 

A  local  modification  to  the  structure  is  defined  as  a  modification  with 
a  p/N  ratio  much  less  than  1.  Appendix  II  demonstrates  that  a  form 
of  the  modified  system's  characteristic  equation  is 


(57) 


g(\  )  =  determinant  (I  +  S(A  )  F(A  ))  =  0 
®  n  -p  -  n  -  n 

where 

zeroes  of  g(A)  (eigenvalues  of  the  modified  structure) 

S(A)  =  S  / j  i  j2  ...  jp\  pxp  (58) 

\jl  h  jp) 

=  condensed  modification  matrix 

A 

F(X)  =  I  /jj  j2  •••  jp\  P*P  (59) 

U  j2  •••  V 

=  condensed  receptance  matrix 

I  =  unit  matrix  of  order  p 
“P 

Equations  (58)  and  (59)  utilize  the  notation  of  equation  (4«).  As  seen 
above,  the  generalized  receptance  matrix  enters  into  the  reanalvsis  pro- 

A 

cedure  through  the  matrix  F(A).  The  zeroes  of  g (A )  are  obtained  iteratively, 
using  either  the  Newton  Raphson  or  Muller  approach.  During  the  iteration 

A 

the  condensed  receptance  matrix  F(\)  must  be  repeatedly  evaluated  for  each 
guess  of  A.  In  addition  a  search  for  an  optimal  structural  design  requires 
computation  of  the  eigenvalues  for  each  trial  design.  These  designs  result 
in  the  formation  of  AH,  AK,  and  AH  modification  matrices  and  consequently 
in  the  formation  of  condensed  effective  modification  matrices  S_.  This 

A 

emphasizes  the  requirement  for  computation  of  F(A)  for  many  values  of  A. 

The  eigensolution  reanalvsis  procedure  is  therefore  a  reasonable  application 
for  the  approximate  receptance  formulas  previously  derived. 

NUMERICAL  EXAMPLE 

The  uniform  beam  model  in  Figure  1  can  be  used  to  demonstrate  the 
accuracy  and  efficiency  of  the  receptance  approach  to  an  eigenvalue  reanalysis 


x;  ' 


-  1  sSiatesaaagk  ■  a  -  ********* 


as  outlined  in  Appendix  II.  This  NASA  model  is  a  crude  continuum  repre¬ 
sentation  of  a  large  space  truss.  The  objective  is  to  determine  the 
optimal  location  for  an  "active"  damper  (c)  to  dampen  transient  vibrations 
resulting  from  maneuvers,  dockings,  and  other  dynamic  loadings. 

Classical  root  locus  plots  are  constructed,  varying  the  damper 
strength  (c)  from  0.1  to  4.0  lb.  sec.  in.^  (0.0175  <_  c  <_  0.70  kN  sec.  m  ^) , 
in  0.1  lb.  sec.  in.^  (0.0175  kN  sec.  m  increments.  For  each  damper 
value  the  lower  portion  of  the  beam's  eigensvstem,  as  determined  from 
the  imaginary  part  of  the  eigenvalues,  is  computed.  These  plots  are 
generated  for  various  damper  locations  and  consequently  may  mandate  the 
computation  of  the  partial  eigensvstem  for  hundreds  of  cases. 

In  this  particular  example  the  system  has  50  dof,  i.e.,  N  =  50 
and  the  generalized  displacement  vector  is 

A  =  (Xx  Q1  X2  02  ...  X25  B25)T 

Suppose  a  damper  is  added  at  node  11,  then  o  =  1  and  J_  =  (j^)  =  (21). 

The  general  equations  of  the  previous  section  reduce  to 

A  A  A  AAA 

SU)  =  A  AM  +  AAC  +  AK  =AAC 


AM  =  AK  =  0,  AC  =  (AC)21  n  =  c 
Also,  for  this  case 

—(A)  “  (F)21,21  =  F21,21 

For  this  case,  equation  (75)  of  Appendix  II  yields  the  following  form 
for  the  modified  system's  characteristic  equation 
g ( A)  *  1  +  Ac  F21  2^(A)  =  0 


Since  the  initial  system  is  undamped  and  is  described  by  symmetric 
M  and  K  matrices,  the  receptance  F^  ^  is  computed,  for  each  guess 
of  A,  by  using  equations  (46)  and  (47) 


R  -  0 


F21,21(A) 


*> 


k=l  m^(.\“+J2^) 


R  =  1 


F21,21(X) 


<^1)21,21 


21, k 


k=i  \cia2«£) 


where  typically  L  <  <  N,  i.e.  an  incomplete  eigensystem  is  used,  and 
$  ,  bk  and  (2.  are  modal  properties  of  the  undamped  system.  The  zeroes 

of  g(A)  are  eigenvalues  of  the  damped  system,  and  are  obtained  with 
Muller's  method.  Root  locus  plots  for  a  damper  at  node  11,  are  shown 
in  Figure  2.  Optimum  damper  and  optimum  damping  ratio  (£)  values  are 
listed  for  each  mode.  Twelve  (L  =  12)  original  system  modes  were 
employed  in  the  receptance  formulas  ,,^),  above.  Figure  2  shows 

results  for  both  the  R  =  0  and  R  =  1  receptance  formulas.  The  eigenvalue 
solver  EISPACK^  was  used  to  judge  the  accuracy  of  the  approximate  re¬ 
analysis  results.  EISPACK  conditions  the  input  matrix  (A^  B,,  see 
equation  (13))  with  subroutine  RAI.ANC,  reduces  the  conditioned  matrix 
to  Hessenberg  forraA^  with  subroutine  ELMHES,  and  then  computes  the 
eigenvalues  of  this  similar  matrix  with  the  QR  algorithm,  utilizing  sub¬ 
routine  HQR.  Table  I  shows  the  maximum  percent  error  for  each  of  the 
eight  (8)  lowest  modes,  for  damper  values  in  the  set  (0.1,  0.2,  0.4, 

0.7,  1.1,  1.2,  1.4,  1.5,  4.0}  lb  sec.  inT1  (0.0175,  0.035,  0.07,  0.1925, 
0.21,  0.245,  0.2625,  0.70  kN  sec.  at  A) .  Note  by  comparison  with  Figure 


2  that  this  set  of  (c)  values  encompasses  all  of  the  optimum  dampers, 
and  the  minimum  and  maximum  damper  values.  The  results  clearly  show 


that  the  computation  of  the  receptances  bv  the  R  *  1  formula  provides 
substantially  improved  eigenvalue  reanalysis  accuracy,  as  compared  to 
the  R  »  0  receptance  computation  results.  Table  II  shows  a  computa¬ 
tion  time  efficiency  comparison  for  the  exact  (F.ISPACK)  and  approximate 
(reanalysis)  methods.  A  large  computer  time  reduction  is  seen  to  result 
with  the  approximate  reanalysis  method.  For  instance,  the  average  re¬ 
analysis  time  per  value  of  C  for  the  R  =  1  approach  is  0.4725  cp.  secs., 
while  the  same  computation  requires  37.6  cp.  secs,  with  EISPACK.  This 
represents  a  computation  time  reduction  by  a  factor  of  79.6. 

This  numerical  example  only  includes  eigenvalue  computations  since 
the  root  locus  plots  do  not  require  eigenvectors. 

CONCLUDING  REMARKS 

A  representation  for  the  approximate  calculation  of  the  receptance 
elements  for  a  general  second  degree,  square  lambda  matrix  has  been 
presented.  The  computational  efficiency  of  solving  for  a  submatrix  of 
the  receptance  is  increased,  in  particular  when  the  submatrix  is  evaluated 
for  many  distinct  values  of  A.  The  subraatrices  are  useful  in  eigensolution 
reanalysis,  design,  or  component  synthesis  schemes,  as  well  as  in  forced 
harmonic  response  analyses.  The  accuracy  of  the  receptances  are  shown 
to  be  improved  relative  to  a  classical  spectral  representation  when  em¬ 
ploying  only  an  incomplete  eigensystem.  Previous  receptance  formulas 
were  shown  to  be  special  cases  of  the  general  results  derived.  These 
prior  works  apply  only  to  undamped,  nonrotating  structures,  with  symmetrical 
property  matrices. 

It  was  assumed  in  equations  (18)  -  (20)  that  the  eigenvalues  of  the 


structural  model  are  all  distinct.  In  addition,  in  equation  (63)  the 


assumption  was  made  that  no  eigenvalue  is  left  invariant  during  a  structural 
modification.  Because  of  the  lack  of  ideal  symmetry  characteristics,  both 
of  these  conditions  are  satisfied  in  many  structural  models.  These 
limitations  on  the  applicability  of  the  formulation  are  subjects  of  ongoing 
research. 


APPENDIX  I :  NOTATION 


General  Notation 

=  matrix-vector  quantity 

A  *  =  inverse  of  matrix  A 

A>0  =  A  symmetric,  positive  definite 

A  >  0  =  A  symmetric,  positive  semi-definite 

diag  (  )  =  diagonal  matrix 

(  )  =  complex  conjugate 

6..  =  Kronecker  delta 

ij 

Re  (  )  =  real  part  of 

i  =  V- 1 
T 

A  =  transpose  of  matrix  A 

Aj,  Bj,  A^,  =  first  order  forms  of  system  property  matrices 

a...  a„.  =  modal  normalization  constants 

It  2i 

C  =  damping  matrix  NXN 
£_(\)  =  a  submatrix  of  _F ,  equation  (48) 

F  =  generalized  receptance  matrix 

G  =  classical  receptance  or  dynamic  flexibility  influence 
coefficient  matrix 
n  =  stiffness  matrix  NXN 

L  =  number  of  system  eigensolutions  available  for  use  in  a 
spectral  type  representation  of  F  or  G 
M  =  mass  matrix  NXN 


modal  normalization  constant 


m 


i 


N  = 
P  = 

Q  = 
R  = 
V  = 


1 

6. 


k 


A.  = 
-1 


= 


total  number  of  system  model  degrees  of  freedom 

2N  X  2N  matrix  pencil 

N  X  2N  right  eigenvector  modal  matrix 

receptance  higher  mode  factor 

N  X  2N  left  eigenvector  modal  matrix 

ith  complex  eigenvalue  of  a  general  vibratory  system 

2N  X  2N  modal  matrix  of  right  eigenvectors  from  a 

first  order  system  formulation 

kth  column  of  P 

2N  x  2N  modal  matrix  of  left  eigenvectors  from  a 
first  order  system  formulation 
kth  column  of  ^ 

ith  modal  damning  constant  for  a  proportionately 
damped  system,  with  symmetric  M,  K,  C 
ith  modal  damping  constant  for  a  symmetric, 
proportionately  damped  system 

ith  complex  right  eigenvector  of  a  general  vibratory 
system 

a  complex  scalar 

nth  generalized  coordinate  vector  for  a  modified 
structural  system 
ith  element  of  £ 


♦u 


j th  eigenvector  of  an  undamped  vibratory  system 
model,  with  symmetric  M,  K 


ith  component  of  <}k 


nth  right  eigenvector  of  a  modified  structural 
sys  tem 

forcing  frequency 

jth  eigenvalue  of  an  undamped  vibratory  system 
model,  with  symmetric  M,  K 


APPENDIX  II:  REANALYSIS  THEORY 


The  free  vibration  problem  for  a  modified  structural  system  is 
expressed  by 

(\2(M  +  AM)  +  An(C  +  AC)  +  (K  +  AK))i^  =0  (N  x  1)  (60) 

An  effective  modification  matrix  mav  then  be  defined  as 


S(\)  =  \"AM  +  AAC  +  AK 


(61) 


Assume  that  matrix  is  is  very  sparse.  This  sparsity  is  quantified  through 
the  integer  set 


-  =  (jl  j2 


••  V 


(p  x  1) 


(62) 


La, 


where  j  ^  are  the  non-null  row  (or  column)  numbers  of  SI.  A  local 
modification  to  the  structure  is  defined  as  one  which  has  a  P/N  ratio 


much  less  than  1.  Define  the  generalized  receptance  matrix  as 
F(A)  =  (A2M  +  AC  +  K)_1 


(63) 


which  exists  as  long  as  A  is  not  an  eigenvalue  of  the  unmodified  system. 
Substitute  the  change  of  basis 

^  =  I(*n)£n(*n'>  (N  x  1)  (64) 


into  equation  (60).  This  yields 


%  *  siyza,))!,  -  2 


By  the  definition  in  equation  (62)  it  follows  that 


(N  x  1) 


-  o. 


k  i  (j,  U  ...  j„) 


and 


*kn  “  -  VJ1  J 2  -p- 


n  +  -j  ^  =  °’  r  e(1’  2’  p) 

j  r  j  r 


■:  -'-tm &**&<!*&*  •, 


(65) 


(66) 


(67) 


1 


where 


S .  =*  i  row  of  S 

— i  J  r  — 

J  r 


(N  x  p) 


r 

^j->n 


'I 

n’ 

P 


T 


(68) 

(69) 


(70) 


Taking  advantage  of  che  sparsity  of  the  columns  of  S_,  equation  (67) 
condenses  further  to 

(Ip  +  i(Xn)F(Xn))Cn  =0  (p  x  1)  (71) 

where  (in  the  notation  of  equation  (48)) 


I-dj/ji  J2 
\jl  j2 

Equation  (71),  accompanied  by 


(p  X  p) 

(72) 

M 

(p  X  p) 

(73) 

j  J 

D/ 

(N  x  1) 

(74) 

represents  the  final  condensed  form  of  the  quadratic  eigenvalue  problem 
(60)  for  the  modified  system.  The  rank  of  the  coefficient  matrix  in 
equation  (71)  must  be  less  than  p,  for  a  nontrivial  null  space  to  exist. 
This  implies  that 


j 

i  I 

i 


i 


g(.\  )  =  dec  (I  +  S(A  )F(A  ))  =  0 
n  — p  —  n  —  n 


(75) 


The  zeroes  of  g(A  )  are  eigenvalues  of  the  modified  structural  system. 

n  J 

Consequently,  g(\)  is  a  form  of  the  modified  system's  characteristic 
equation . 
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Table  I  Accuracy  Comparison:  Approximate  (Reanalysis) 
Versus  Exact  (EISPACK)  Eigenvalues  for  a  Damped 
Beam* 


R  =*  0 

Re  (A) 

Max  %  error 

R  =  1 

Re  (A) 

Max  %  error 

R  =  0 
lm(A ) 

Max  %  error 

R  =  1 

Im(A ) 

Max  /  error 

X1 

-0.008 

0.000 

0.000 

0.000 

X2 

-0.016 

0.000 

0.000 

0.000 

X3 

-0.073 

0.000 

0.016 

0.000 

X4 

-0.458 

0.000 

0.042 

0.000 

X5 

-1.440 

-0.005 

0.125 

0.000 

X6 

-2.17 

-0.024 

0.248 

0.003 

X7 

-9.83 

-0.318 

0.016 

0.000 

-7.82 

-0.387 

0.804 

0.038 

*  r 
Data  shown  is  the  maximum  percent  error  with  C  e  10. 1,  0.2,  0.4,  0.7, 

1.1,  1.2,  1.4,  1.5,  4.0}  lb.  sec.  inT1  (0.0175,  0.035,  0.07,  0.1925, 

0.21,  0.245,  0.2625,  0.70)  kN  ?ec.  m'1 


i 


* 


r  j 
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Table  II  Computation  Time  Comparison:  Approximate  (Reanalysis) 
Versus  Exact  (EISPACK)  Approaches 


(a)  Formation  of  K,  M  and  computation  of  original  system  eigensolutions : 
19.46  cp.  secs. 

(b)  R  =  0  reanalysis  for  c  e  {0.1,  0.2,  ...  4.0},  40  values  of  c  and  10 
A^  per  value:  16 . 8  cp.  secs.,  average:  0.42  sec/damper  value 

(c)  R  =  1  reanalysis  for  c  c  {0.1,  0.2,  ...,  4.0},  40  values  of  c  and  10 
A^  per  value:  18.9  cp.  secs.,  average:  0.4725  sec/danper  value 

(d)  EISPACK  solution  for  9  values  of  c  and  10  A.  per  value  c  £  {0.1,  0.2, 
0.4,  0.7,  1.1,  1.2,  1.4,  1.5,  4.0}:  338 . 0  cp.  sec.,  average:  37.5556 
sec/damper  value 


K  =  0.10  lh/in  (17.51  N/m) 

H.  =  6.0  in.  (15.24  cm) 

J 

T.  =  0.1375  in.  (0.47625  cm) 

J 

L.  =  6.0  in.  (15.24  cm) 

L  =  144  in.  (3.6576  m) 

pg  =  0.0967  lb/in3  (0.02596  N'/cm3) 

E  =  10  x  106  lb/in2  (6S.93  GN/m2),  Ailuminum 

0.1  <  c  <  4.0  lb. sec. /in  (0.0175  <-  c  <  0.70  I'.N. sec/m) 


Original  system  natural  frequencies:  0^  =  2.065  rad/sec,  =  3.804, 
V.3  =  12.504,  =  32.472,  9.^  =  63.211,  =  104.332,  9^  =  155.792, 

Q8  =  217.587 


Figure  1  Euler  Bernoulli  beam  model  for  j 

space  structure  reanalvsis  1 
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MODE 

OPTIMUM  DAMPER 

OPTIMUM  ZETA 

1 

.100 

.71313 

2 

.100 

.01470 

3 

.200 

.18230 

4 

.400 

.07575 

5 

.700 

. 07624 

£ 

.700 

.10437 

7 

1.400 

.00130 

8 

1.100 

.08778 

MODE 

OPTIMUM  DAMPER 

OPTIMUM  ZETA 

1 

.100 

.71322 

2 

.100 

.01470 

3 

.200 

. 18243 

4 

.400 

.07602 

5 

.700 

.07723 

6 

.700 

. 10665 

7 

1.500 

.00201 

8 

1.200 

.03508 

RUN  NO.  40781.200 
DAMPER  NOOK  NO.  11 

NO.  DAMPER  VALUES  40  . 

SMALLEST  DAMPER(*)  .100  LB . SF.C .  IN_1  (0.0175  kN.SEC.m~1) 

DAMPER  INCREMENT  .100  LB.SEC.IN~1  (0.0175  kN.SEC.m~1) 
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Figure  2  Root  locus  plot  for  space  structure  beam 
model  utilizing  eigenvalue  reanalysis 


